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We present the first measurements of cfiarge-dependent correlations on angular difference vari- 
ables ?7i — ?72 (pseudorapidity) and (j)i — (j)2 (azimuth) for primary charged hadrons with transverse 
momentum 0.15 < Pt < 2 GeV/c and jr;| < 1.3 from Au-Au collisions at ^ysNN ~ 130 GeV. We ob- 
serve correlation structures not predicted by theory but consistent with evolution of hadron emission 
geometry with increasing centrality from one-dimensional fragmentation of color strings along the 
beam direction to an at least two-dimensional hadronization geometry along the beam and azimuth 
directions of a hadron-opaque bulk medium. 

PACS numbers: 24.60.Ky, 25.75.Gz 



I. INTRODUCTION 



II. ANALYSIS METHOD 



Analysis of correlations and fluctuations plays an im- 
portant role in studies of the colored medium produced in 
ultrarelativistic heavy-ion collisions jl], |^ y| . In-medium 
m,odification of parton scattering and fragmentation of 
energetic partons by the bulk medium produced in heavy- 
ion collisions may significantly alter large-momentum- 
scale two-particle correlations relative to those observed 
in p-p collisions. Large-momentum-scale correlations 
may result from initial-state multiple scattering [^ 
in-medium dissipation of scattered energetic partons 
and hadronization of the colored medium to final-state 
hadrons (fragmentation of color strings in p-p, hadroniza- 
tion of the bulk medium in A- A). The local geometry of 
hadronization, which can be accessed by net-charge cor- 
relations, is the subject of this paper. 

String fragmentation models |3| describe two-particle 
correlations on pseudorapidity and azimuth (77, <f)) in 
high-energy p-p collisions in terms of local conservation of 
transverse momentum and net charge leading to canon- 
ical suppression of event-wise net-momentum and net- 
charge fiuctuations. The nature of the corresponding pro- 
cess in A-A collisions remains an open question. Some 
change should be expected in the correlation structure 
as the medium evolves from that produced in very pe- 
ripheral collisions (approximating minimum-bias proton- 
proton collisions) to that in central heavy-ion collisions. 
Predictions have been made of dramatic suppression of 
net-charge fluctuations in central A-A collisions as a 
signal of quark-gluon plasma formation 8J. The ques- 
tion arises what detailed net-charge correlation structure 
would correspond to such predictions, and what structure 
is actually present in heavy-ion collisions. 



In this Letter we report the first measurements in 
heavy-ion collisions of the centrality dependence of two- 
particle charge- dependent (net-charge) correlations on 
angular subspace (ry, 0), where charge-dependent here 
refers to the difference between correlations for like- 
charge-sign pairs and unlike-sign pairs. This analysis 
is based on Au-Au colhsions at y/SNN = 130 GeV ob- 
tained with the STAR detector at the Relativistic Heavy 
Ion Collider (RHIC). The observed correlation structure 
suggests that local charge conservation at hadroniza- 
tion combined with increasing system density and spatial 
extent results in evolution with Au-Au centrality from 
one-dimensional (ID) charge- ordering (locally alternat- 
ing charge signs) on configuration space z (the collision 
axis), coupled to pz (or pseudorapidity vj) by longitudi- 
nal Bjorken expansion, to two-dimensional (2D) charge 
ordering on beam and azimuth directions (z,;/)). Those 
results have not been anticipated by theoretical mod- 
els nn. 



We wish to access the complete charge- dependent (CD) 
structure of two-particle density p{pi,p2) with minimal 
distortion and without imposition of a correlation model. 
In this analysis of net-charge angular correlations we 
project the two-particle momentum space onto angular 
subspace (?yi,?/2, 0i,</'2) by integrating over a specific 
transverse momentum interval. The structure of net- 
charge correlations on transverse momentum with spe- 
cific angular constraints will be considered in a future 
analysis. 

Correlations are obtained with a differential analysis 
which compares object and reference pair density distri- 
butions. The object distribution is comprised of parti- 
cle pairs formed from single events, referred to as sibling 
pairs, and the reference distribution consists of pairs com- 
bining particles from two different but similar events, re- 
ferred to as mixed pairs. The corresponding pair densities 
are denoted by Psib{pi,P2) and Pmix{pi,P2) respectively. 
The two-particle correlation function C (as commonly 
defined in nuclear physics) and pair-number density ra- 
tio r (as used in the study of quantum correlations or 
HBT [13) are then defined and related by 



C (pi , P2 ) = Psib{Pl,P2)~ Pmix {Pl , P2 ) 
= Pinix{Pl,P2){r{pi,P2)-'^), 



(1) 



with r = Psib/Pmix- Difference r — 1 is the correlation 
measure we use. In order to visualize the CD correla- 
tion structure in the 4D angular subspace (?7i, ?72, 4>i,4'2) 
pair densities can be projected onto separate 2D sub- 
spaces (771,772) and {4>i,4>2)- Those projections, discussed 
further below, discard a substantial amount of the in- 
formation in the full two-particle space. However, they 
reveal that significant variation is restricted to difference 
variables r^^ = rji — rj2 and 0a = 4>i ^ 4>2 (the notation 
is explained in Sec. lIVp . For this analysis we therefore 
simultaneously project the 4D subspace onto those angu- 
lar difference variables. The resulting 2D distribution is 
referred to as a joint autocorrelation. An autocorrelation 
is a projection by averaging 11] from subspace (a;i,a;2) 
onto difference variable xa = xi — a;2. A joint autocor- 
relation is a simultaneous projection onto two difference 
variables. The result of this projection technique is a 
nearly lossless (distortion free) projection from the initial 
4D angular subspace onto a 2D autocorrelation space. 

In this analysis, sibling and mixed pair- number 
densities p{pi,p2) for four charge-pair combinations 

(-f-l-, -I — , — h, ) were projected onto (771,772), {(j)i,(j)2) 

and (?7a,0a)- The projection was done by fill- 
ing histograms of pair numbers riab — ^x^y pi^aiVb), 
where subscripts ab denote the 2D bin indices and 
CxT^y are histogram bin widths on variables x^y S 
{771, 7^2, 01, 02, ??A, 0a}- Sibling and mixed pair- number 
histograms for each charge-pair combination were sepa- 
rately normalized to the total number of detected pairs 



in each centrality class: 



nab,: 



^ab.sib 



I A^ab 



riab,: 



and nab,mix = nab,mix/J2ab'^o.b,mix- Normalized pair- 
number ratios -Tab = f>'ab.siblf''ab,mix ai'C the basis for this 
analysis. 

To reduce systematic error, ratio histograms were ob- 
tained for subsets of events within a given centrality 
class which have similar multiplicities (differences < 50) 
and primary collision vertex locations within the detec- 
tor (within 7.5 cm along the beam axis). Ratios fab for 
each centrality class were defined as weighted (by total 
number of sibling pairs) averages over all subsets in that 
centrality class. Ratios were further combined to form 

like-sign (LS: ++, ), unhke-sign (US: H — ,— +), and 

charge-dependent (CD = LS — US) ratios. In this analy- 
sis we adopt a CD sign convention compatible with stan- 
dard particle physics isospin convention and net-charge 
fluctuation measures |13 |. 



III. DATA 

Data for this analysis were obtained with the STAR de- 
tector 13'] using a 0.25 T uniform magnetic field parallel 
to the beam axis. A minimum-bias event sample required 
coincidence of two Zero-Degree Calorimeters (ZDC); a 0- 
15% of total cross section event sample was defined by 
a threshold on the Central Trigger Barrel (CTB), with 
ZDC coincidence. Event triggering and charged-particle 
measurements with the Time Projection Chamber (TPC) 
are described in [13. Tracking efficiencies, event and 
track quality cuts and primary-particle definition are de- 
scribed in [13, llJI . Charged particles were accepted in 
1 77! < 1.3, full azimuth and transverse momentum [pt) 
range 0.15 < pt < 2 GeV/c. Particle identification was 
not implemented but charge sign was determined. Cor- 
rections were made to ratio f for two-track inefficiencies 
due to overlapping space points in the TPC (merging) 
and intersecting trajectories reconstructed as > 2 parti- 
cles (splitting) by applying two-track proximity cuts in 
the TPC to both psib and Pmix similar to that done in 
HBT analyses. 

Small-momentum-scale correlation structures due to 
quantum. Coulomb and strong-interaction correla- 
tions [lOj were suppressed by eliminating sibling and 
mixed track pairs (~22% of total) with |77a| < 1.0, 
|(/ia| < 1-0 and \pti-pt2\ < 0.2 GeV/c iipt < 0.8 GeV/c 
for either particle. Those cuts do not significantly af- 
fect the correlation structures shown here. Four central- 
ity classes for 300k events labeled (a) - (d) for central 
to peripheral were defined by cuts on TPC track multi- 
plicity N within the acceptance defined here relative to 
minimum-bias event multiplicity frequency distribution 
upper half-maximum end-point Nq, which corresponds 
to the maximum participant number |l2l |l5J . Four cen- 
trality classes were defined by (d) 0.03 < iV/7Vo < 0.21, 
(c) 0.21 < N/No < 0.56, (b) 0.56 < N/Nq < 0.79 and 
(a) N/No > 0.79. 




FIG. 1: Normalized LS pair-number ratios f for collisions in 
centrality class (a) (most central) for (771,772) (left panel) and 
{(t>i,<l>2) (right panel). 



IV. TWO-PARTICLE DISTRIBUTIONS 

Fig. n shows ratio histograms fab for the LS charge 
combination on (771,772) and (0i,02) for the most cen- 
tral event class, denoted (a). Deviations from unity 
(f — 1) of this per-pair correlation measure contain a di- 
lution factor [l6l | 1/iV {N is defined as the mean mul- 
tiplicity in the detector acceptance) and are therefore 
numerically a few permil for central Au-Au collisions. 
However, the correlation structure is large compared to 
statistical errors (c/. Figs. ElUl . A sinusoid associated 
with elliptic flow (consistent with conventional reaction- 
plane measurements) dominates the ((/)i, 02) correlations 
in the right panel. The a7if? correlated LS distribution 
on (771,772) in the left panel (anticorrelated: depression 
along the 771 = 772 diagonal) suggests charge ordering 
from longitudinal string fragmentation as in p-p colli- 
sions iLdZiJ- However, these correlations projected sepa- 
rately onto (771, 772) and {(j)i,(j)2) are incomplete, and quite 
misleading for A- A collisions. A more complete picture 
is obtained from 2D joint autocorrelations on difference 
variables (r]/^,(J3^) as shown in Fig. |21 

Because of the symmetry of these distributions on the 
angular spaces (a;i,a;2) their description is more natural 
on diagonal sum and difference variables xy: and xa (re- 
serving conventional difference notation Aa; for displace- 
ment on a ID space x). The invariance of correlation 
structure on sum variables rjs = ?7i +^72 and c/jy: = (j>i+<j)2 
in Fig. n (i.e., parallel to the 771 = 772 or 4>i = 4>2 diago- 
nals) implies that each distribution can be projected onto 
its difference variable 0a = 01 — 02 and rj^ = r]i — 772 
to form an autocorrelation without loss of information. 
The projection is done by averaging bin contents along 
each diagonal in Fig.nparallel to the sum axis {e.g. the 
771 = 772 diagonal) to obtain the bin contents of a ID 
autocorrelation on 77A or 0a (the difference axes). Auto- 
correlation details are described in |l§ . ll9l |. If projections 
are made simultaneously onto both difference variables of 
Fig. nthe resulting 2D joint autocorrelation on (77A, 0a) 
compactly represents all signiflcant correlation structure 
on 4D angular subspace (771,772,01,02). 

In Fig. 121 perspective views are shown of CD joint au- 
tocorrelations for four centrality classes of Au-Au col- 
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FIG. 2: Perspective views of two-particle CD joint autocor- 
relations iV(f — 1) on {rjA, 4>a) for central (a) to peripheral (d) 
collisions. Center bins at (j!)a = f?A = 0, containing photon- 
conversion electron pairs, were omitted from model fits. 




FIG. 3: Projections of 2D CD autocorrelations N{r — 1) in 
Fig.|2|onto individual difference variables tja (solid dots) and 
(l>A (open triangles) for central (a) to peripheral (d) collisions. 
Solid (dashed) curves represent projections of 2D analytical 
model fits to data on t^a ((/"a). The 2D negative peaks are 
substantially reduced in amplitude after projecting onto ID. 



ERRORS 



lisions at y/s^N = 130 GeV. Quantity N{f - 1) |20| 
represents per-particle correlations {i.e., distribution of 
average numbers of correlated pairs per final-state par- 
ticle) and is 0(1) for all centralities. Distributions in 
Fig. [21 are dominated by a 2D negative peak which is 
broader and elliptical for peripheral collisions (d) with 
major axis along (/)a, transitioning smoothly to a nar- 
rower and deeper peak symmetric on (?7a,0a) for cen- 
tral collisions (a). The negative peak means that unlike- 
sign charge pairs are more probable than like-sign pairs 
for small angular separations on pseudorapidity and az- 
imuth, consistent with local charge conservation (sup- 
pression of net-charge fluctuations). The vertical axis 
limits common to all panels were chosen to enhance the 
visibility of structure at large angular separations as op- 
posed to showing the full depth of the negative peak at 
0A = '7A = 0. Note that no CD (charge-dependent) 
component of elliptic flow is observed at the sensitivity 
level of these data. ID projections of Fig. [21 distributions 
and their 2D model fits (discussed below) onto individ- 
ual difference variables (/)a and tja are shown in Fig. [31 
Solid dots and curves (open triangles and dashed curves) 
correspond to rjA i4>A) projections. The projections are 
over the pair acceptances apparent in Fig. [21 



Statistical errors for f in Fig. ^ (central collisions) are 
±0.00015 for all bins. Statistical errors for ID autocorre- 
lations arc uniform on (j)/^ (since </> is a periodic variable) 
but approximately double as |77a| increases from to 2 
(due to finite 77 acceptance). Statistical errors at t^a ~ 
vary from ±0.00015 for central collisions to ±0.0007 for 
peripheral collisions, again reflecting the 1/N dilution 
factor. In contrast, statistical errors for N{f—1) in Fig. [21 
are approximately ±0.2 (one tick) for 77A ^ and are in- 
dependent of centrality Statistical errors for projections 
in Fig.[31are shown explicitly in that figure by error bars. 
Systematic errors were estimated as in 12j. System- 
atic uncertainties associated with two-track inefficiency 
corrections and small momentum scale correlation cuts 
are negligible for this analysis. Systematic error due to 
non-primary backgrounds (dominant source) 14], whose 
correlation with true primary particles is unknown, is 
estimated to be at most ±7%, assumed uniform for all 
(??a,<^a) in the STAR acceptance. Contributions from 
resonance (p°,tf) decays are estimated to be at most 
about 10% of the negative peaks at (Jja = ?7a = in 
Fig.[21in the range |?7a| < 0.5, |0a| < 2^. 



VI. MODEL FITS 

The distributions in Fig. 13 and their counterpart for 
p-p colHsions [23| reveal two asymptotic forms at the cen- 
trahty Hmits: a ID gaussian on rjA {uniform on 0a) for 
p-p coUisions and a 2D exponential on {rjA, ^a) for cen- 
tral Au-Au collisions. The two forms may be limiting 
cases of a single evolving structure, or they may corre- 
spond to two independent correlation mechanisms with 
complementary centrality trends. A preliminary fitting 
exercise indicated that these 130 GeV Au-Au data do not 
have sufficient statistical power or centrality range to ex- 
plore the possibility of a single evolving peak structure. 
We therefore used the simpler superposition model. 

The distributions in Fig. |21 were fitted with a five- 
parameter model function consisting of a 2D exponential 
function peaked on both ?7a and 4>a and a ID gaussian 
on ?7a, constant on d)/\ (the latter motivated by the p-p 
limiting case [l3;l22l) plus a constant offset, all defined 
relative to quantity f — 1 as 



F = Aq + Ai exp 



A2 exp ■ 




2 / s 2 
A \ , / VA 



(2) 



F interpolates between the ID gaussian peak observed 
in p-p and the 2D exponential peak observed in central 
Au-Au collisions. Correlations between amplitudes Ai 
and A2 were negligible because of the distinct one- and 
two-dimensional peak shapes. Parameters a^^ and cr^^ 
are the r.m.s. widths of the 2D exponential peak when 
projected onto the respective difference variables. 

Best- fit values for varied parameters and x^/DoF for 
the four centralities are listed in TableQ] The width of the 
ID gaussian, most evident near |(/)a| ~ tt in Fig. 2(d), was 
best determined by those peripheral data to be 1.5 ±0.25 
and was held fixed at that value for the other centralities 
to obtain the amplitude estimates. The observed periph- 
eral Au-Au 0A width is definitely larger than the corre- 
sponding width for p-p collisions. Also included is track- 
ing efficiency-correction factor S [23. Total systematic 
error for efficiency-corrected amplitudes in Table ^ was 
11% (errors added in quadrature). The model fits indi- 
cate that with increasing centrality the 2D exponential 
peak exhibits 1) strong amplitude increase, 2) significant 
width reduction and 3) approach to approximately equal 
widths on (f>A and tja for central collisions (c/. Fig. 13 
e.g., at mid-rapidity a,-i^ — 0.6 corresponds to polar an- 
gle difference 0.57, which is directly comparable to cr^^). 



VII. DISCUSSION 



TABLE I: Parameters and fitting errors (only) for model fits 
[Eq. ^] to joint autocorrelation data in Fig.|21for centrality 
bins (a) - (d) (central - peripheral). Total systematic error 
for tracking efficiency-corrected amplitudes is 11% ^3 



centrality 


(d) 


(c) 


(b) 


(a) 


error'"(%) 


S [23] 

N 


1.19 
115.5 


1.22 
424.9 


1.25 
789.3 


1.27 
983.0 


8 (syst.) 


SNAo 


0.98 


0.80 


0.91 


0.79 


11-12 


SNAi 

O"0A 


-4.1 
0.94 
0.66 


-6.8 
0.75 
0.59 


-7.7 
0.72 
0.58 


-7.7 
0.72 
0.58 


6-4 
11-5 
10-5 


SNA2 


-0.51 


-0.11 


-0.15 


-0.021 


0.17-0.19 * 


xVdof 


380 
315 


315 
3157 


314 
31 5 


329 
315 





"Range of fitting errors in percent, from peripheral to central. 
''Magnitude of fitting errors. 



Au collisions differ dramatically from those for p-p col- 
lisions. CD angular correlations for p-p collisions are 
dominate d by a ID negative gaussian peak on tja with 
^»?A — 1 [l3i ^M ' conventionally associated with longi- 
tudinal charge ordering on z during string fragmenta- 
tion 7ji plus a 2D gaussian peak associated with quan- 
tum correlations. For the most peripheral Au-Au cen- 
trality (d) in this analysis we observe CD correlation 
structure intermediate between p-p and central Au-Au 
collisions, consistent with the fact that collision events 
in centrality class (d) for these 130 GeV data are not 
very peripheral: they contain about 100 particles in the 
STAR acceptance (see Table IJ). In central Au-Au colli- 
sions the ID gaussian peak is no longer detectable. In- 
stead, a large-amplitude 2D negative exponential peak 
dominates the correlation structure, with similar widths 
on ?7A and (f)A much reduced from those measured in p-p 
collisions. 
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FIG. 4; Left panel: Efficiency corrected correlation ampli- 
tudes for 2D exponential (dots) and ID gaussian (triangles) 
components from Tablelllfor neg ative peaks in Fig.|5]are plot- 
ted on mean path length u [13. Right panel: Fitted widths 
a^^ (dots) and tan~^ a^^ (triangles) are plotted on v. Plot- 
ting variable tan~^ permits the divergent p-p a^^ value to be 
included. Hatched regions and u — 1 data points summarize 
p-p limiting values. Curves guide the eye. 



This analysis demonstrates for the first time that 
charge-dependent angular correlations for central Au- 



Variations of peak amplitudes and widths with Au-Au 
centrality are shown in Fig. ^ along with p-p limiting 



cases (cross-hatched bands) from STAR p-p data at 200 
GeV Q^, consistent with ISR p-p data at 52.5 GeV [l3|. 
The p-p data points in Fig. ^ (values at j^ = 1) indicate 
the aniphtude and r.m.s. width of the ID gaussian on 
r^A, the uniformity of that correlation on c/)/^ (c^^ 3> 1) 
and the absence of a 2D exponential on (77A, <J>a) in the fit 
residuals, represented by the solid dot in the left panel at 
1/ = 1. Comparison of the low-pt (0.15 < pt < 0.5 GeV/c) 
p-p results with the present Au-Au results is qualitative 
but reasonable given the similarity in shape of the Au-Au 
CD correlations for 0.15 < pt < 0.5 (discussed below) to 
those in Fig. |21 

The collision centrality is represented by mean partici- 
pant path length 1/ 15], defined as the average number of 
nucleons encountered by a participant nucleon. That cen- 
trality measure is desirable because it permits compar- 
isons with p-A collisions, initial-state scattering should 
follow a trend linear in 1/ and v also provides an estimate 
(proportionality) of final-state pathlength. 

We adopt the strategy of plotting tan^^ (cr^^) rather 
than (T0^ so as to include the p-p 'infinite azimuth width' 
on the same plot, since that distribution is approxi- 
mately uniform on (j). Interpolations among the measured 
Au-Au points are sketched by the solid and dash-dot 
curves. Extrapolations to corresponding p-p values are 
sketched by the dashed and dotted curves. The extrapo- 
lations contain substantial uncertainties in relating p-p to 
mid-peripheral Au-Au results. Efficiency- corrected per- 
particle correlation amplitudes SNA for central Au-Au 
collisions exceed in magnitude those for p-p collisions by 
a factor 10. The dramatic shape and amplitude changes 
strongly contradict a p-p linear superposition hypothe- 
sis ,2QJ for all but the most peripheral Au-Au collisions. 

These results for net-charge angular correlations sug- 
gest that CD correlations in Au-Au collisions, as in p-p 
collisions, derive from configuration-space charge order- 
ing as a consequence of local charge conservation during 
hadronization, but the hadronization geometry changes 
from ID (ji) in p-p collisions to at least 2D (77, (f>) in cen- 
tral Au-Au collisions, leading to an approach to angular 
symmetry on {ri^,(j)/^). Transverse charge ordering (on 
Pt) is also possible but is studied in a separate analy- 
sis. Hadronic rescattering in A-A collisions could reduce 
the CD correlation amplitude at large 0a but would also 
reduce the width on t^a and therefore cannot be solely 
responsible for the nearly symmetric peak shape in cen- 
tral Au-Au collisions. In Fig.Qlthe contribution from ID 
charge ordering (gaussian peak on j]a) is already sub- 
stantially reduced for centrality (d) {v ^ 2.5) in favor of 
the symmetric component (exponential peak). 

A hadron-opaque medium in more central collisions 
may contribute to the newly-observed exponential peak 
shape. An exponential distribution on pair opening an- 
gle [radius on (77, (j>)] is consistent with: 1) correlations 
detected only if both members of a correlated pair are 
not significantly scattered, 2) scattering probability de- 
termined by a mean free path, 3) mean path length in 
the medium increasing monotonically with pair opening 



angle. That rescattering picture assumes that CD cor- 
relations do not result from hadronization outside the 
medium. Contributions from charge ordering in jet frag- 
mentation were studied by splitting central Au-Au data 
at Pt = 0.5 GeV/c, below which jet fragments should 
be negligible. Negative peak structures as in Fig. [3 were 
observed to dominate both subsamples, although the am- 
plitudes were not identical. 

HIJING and RQMD 0, charge-dependent angu- 
lar correlations qualitatively disagree with data. Hi- 
JING charge-dependent correlations are determined by 
the Lund model 0] via PYTHIA |2J| , and are consequently 
consistent with p-p ID string fragmentation for all A-A 
centralities: a ID gaussian on 77A with amplitude about 
10% of the exponential peak in Fig.|21(a). RQMD, rep- 
resenting mainly resonance decays and hadronic rescat- 
tering, exhibits a broad 2D gaussian on (r7A,</'A)j with 
amplitude also about 10% of the exponential peak in the 
data for central collisions. Large-scale correlations as in 
Fig. ^observed for US and LS pairs in data are consistent 
with local charge ordering but inconsistent with CD cor- 
relations from decays of hadronic resonances such as the 
p*^, which affect only the US pair type. That observation 
further argues against a resonance-gas scenario. 

Measurements of net-charge fluctuations have been ad- 
vocated as a probe of heavy-ion collisions. Predictions 
of dramatic suppression of net-charge fluctuations in the 
case of QGP formation based on entropy arguments [g 
refer by implication to an integral of net-charge angu- 
lar correlations over a detector acceptance. Phenix ob- 
served net-charge fluctuations in Au-Au at 130 GeV [23 
slightly reduced from 'stochastic behavior' and indepen- 
dent of collision centrality. The data were consistent 
with RQMD representing a resonance gas. STAR ob- 
served net-charge fluctuations in Au-Au at 200 GeV z6\ 
intermediate between what is expected from canonical 
suppression in a partial acceptance and a resonance gas, 
again with little or no centrality dependence. Those con- 
clusions are in sharp contrast to what we observe in the 
present analysis. 

It is important to note that net-charge fluctuations 
within a given detector acceptance integrate CD joint 
autocorrelations such as those presented in this paper 
(within a constant offset) over that acceptance, as de- 
scribed in 18] . As integral quantities, fluctuation mea- 
surements are insensitive to the differential structure of 
angular correlations. In the present analysis we observe 
dramatic changes in differential structure (10-fold am- 
plitude increase, nearly two-fold width reduction) while 
corresponding peak integrals exhibit only modest change 
with collision centrality (integrals of observed CD peaks 
using peak parameters in TableQ]increase linearly in mag- 
nitude on ly by about 20%). We suggest that the the- 
oretical connection between net-charge fluctuation sup- 
pression and QGP formation, currently based only on 
large-scale integral measures, should be re-examined in 
the more differential context of CD autocorrelation struc- 
ture. 
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VIII. SUMMARY 

In summary, we have measured charge-dependent an- 
gular correlations on pseudorapidity and azimuth differ- 
ence variables (r/i — 1J2) and {4>i — (j)2) for Au-Au colli- 
sions at y/sNN = 130 GeV. The data are consistent with 
local charge conservation or canonical suppression of net- 
charge fluctuations, evolving from ID (along rj) color- 
string fragmentation in p-p collisions to exponentially- 
attenuated (on opening angle) 2D charge-ordered emis- 
sion from a hadron-opaque medium in central Au-Au col- 
lisions. The transition from ID to 2D correlation struc- 
ture occurs rapidly with increasing collision centrality. 
These results are qualitatively inconsistent with predic- 
tions from standard Monte Carlo collision models typi- 
cally applied to single-particle differential distributions 
and integrated yields from relativistic heavy-ion colli- 
sions. Charge-dependent angular autocorrelations pro- 



vide unique differential access to the changing geometry 
of hadronization and hadronic rescattering as the energy 
density and spatial extent of A-A collisions vary with 
centrality. 

We thank the RHIC Operations Group and RCF at 
BNL, and the NERSC Center at LBNL for their sup- 
port. This work was supported in part by the HENP 
Divisions of the OfSce of Science of the U.S. DOE; the 
U.S. NSF; the BMBF of Germany; IN2P3, RA, RPL, 
and EMN of France; EPSRC of the United Kingdom; 
FAPESP of Brazil; the Russian Ministry of Science and 
Technology; the Ministry of Education and the NNSFC 
of China; IRP and GA of the Czech Republic, FOM of 
the Netherlands, DAE, DST, and CSIR of the Govern- 
ment of India; Swiss NSF; the Polish State Committee 
for Scientific Research; STAA of Slovakia, and the Korea 
Sci. & Eng. Foundation. 



[1] R. Stock, Nucl. Phys. A661, 282c (1999); H. Heiselberg, 

Phys. Rep. 351, 161 (2001). 
[2] A. Dumitru, R. Pisarski, Phys. Lett. B504, 282 (2001). 
[3] L. M. Bettencourt, K. Rajagopal and J. V. Steele, Nucl. 

Phys. A693, 825 (2001). 
[4] M. Gazdzicki, A. Leonidov, G. Roland, Eur. Phys. J. C6, 

365 (1999). 
[5] X.-N. Wang, M. Gyulassy, Phys. Rev. D 44, 3501 (1991). [17] 

[6] Q. J. Liu and T. A Trainor, Phys. Lett. B567, 184 (2003). 
[7] B. Andersson, G. Gustafson, G. Ingelman and T. [18] 

Sjostrand, Phys. Rep. 97, 31 (1983). 
[8] M. Asakawa, U. Heinz, B. Miiller, Phys. Rev. Lett. 85, [19] 

2072 (2000); S. Jeon, V. Koch, Phys. Rev. Lett. 85, 2076 
(2000). 
[9] H. Sorge, H. Stocker, W. Greiner, Nucl. Phys. A498, 

567c (1989); Ann. Phys. (N.Y.) 192, 266 (1989). [20] 

[10] C. Adler et al. Phys. Rev. Lett. 87, 082301 (2001). 
[11] Averaging rather than simple projection is an essential 
feature of autocorrelations required to properly account 
for acceptance effects in two-particle spaces. [21] 

[12] J. Adams et al. (STAR Collaboration), Phys. Rev. C 71, [22] 

064906 (2005). 
[13] K. H. Ackermann et al, Nucl. lustrum. Meth. A 499, [23] 

624 (2003); see other STAR papers in volume A499. 
[14] C. Adler et al, Phys. Rev. Lett. 87, 112303 (2001); ibid. 

89, 202301 (2002). 
[15] Centrality measure u estimates the mean participant [24] 
path length as a number of encountered nucleons. For this 

analysis U = 5.5 {N/NqY^^ ~ 5.5 {Npart/Npart,maxy^''^ ^ 

'iNbin/Npart, bascd on Glauber-model simulations. Npart [25] 
is the number of participants, N^in the number of binary 
collisions, and No is the upper half-maximum endpoint of [26] 
the minimum-bias data distribution plotted as da/dN^'^. 
[16] Under composition of independent (uncorrelated) sub- 



systems {e.g., hypothetical independent N-N collisions 
within A-A collisions) correlations per particle remain 
constant, whereas correlations per pair go as l/rich- An 
exception to that behavior is quantum correlations for 
bosons (pions) where all pairs with momentum differ- 
ence q below some value are correlated, in which case 
r— 1— »lasg— >0, independent of system size. 
D. Drijard et al. (ACCDHW Collaboration), Nucl. Phys. 
B 166, 233 (1980). 

T. A. Trainor, R. J. Porter and D . J. Prindle, J. Phys. G 
31, 809 (2005) fhep-ph/0410182V 

D. J. Prindle and T. A. Trainor, hep-ph/0506173 To ap- 
pear in: Proceedings of the MIT Workshop on Correla- 
tions and Fluctuations in Relativistic Nuclear Collisions, 
Cambridge, Massachusetts, 21-23 April, 2005. 
N{r — 1), measuring correlations per final-state particle 
(typically 0(1) for all centralities) , is invariant with cen- 
trality if A-A collisions are linear superpositions of p-p 
collisions. 

R. Ray and R. Longacre, nucI-ex/0008009 
R.J. Porter and T. A. Trainor (STAR Collaboration), 
hep-ph/0406330 

Extrapolation factors S for 7V(f — 1) provide corrections 
to amplitudes Ao, Ai and A2 for background contamina- 
tion and tracking inefficiency [ijj- Systematic error in S 
was estimated to be ±8%. 

T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994); 
T. Sjostrand, L. Lonnblad, S. Mrenna and P. Skands, 
pp-ph/0308153 

K. Adcox et al. (PHENIX Collaboration), Phys. Rev. 
Lett. 89, 082301 (2002) (nucbex/0203014). 
J. Adams et al. (STAR Collaboration), Phys. Rev. C 68, 
044905 (2003) ( ,nucl-ex/0307007j . 



